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SUMMARY 

Characteristics of the following six flight paths for decelerating 
from a supercircular speed are developed in closed form: constant angle 

of attack, constant net acceleration, constant altitude, constant free- 
stream Reynolds number, and "modulated roll." The vehicles were required 
to remain in or near the atmosphere, and to stay within the aerodynamic 
capabilities of a vehicle with a maximum lift-drag ratio of 1.0 and 
within a maximum net acceleration G of 10 g’s. The local Reynolds 
number for all the flight paths for a vehicle with a gross weight of 
10,000 pounds and a 60° swept wing was found to be about 0.7X10®. 

With the assumption of a laminar boundary layer, the heating of the 
vehicle is studied as a function of type of flight path, initial G load, 
and initial velocity. The following heating parameters were considered: 
the distribution of the heating rate over the vehicle, the distribution 
of the heat per square foot over the vehicle, and the total heat input 
to the vehicle. The constant G load path at limiting G was found to 
give the lowest total heat input for a given initial velocity. For a 
vehicle with a maximum lift-drag ratio of 1.0 and a flight path with a 
maximum G of 10 g's, entry velocities of twice circular appear thermo- 
dynamically feasible, and entries at velocities of 2.8 times circular 
are aerodynamic ally possible. The predominant heating (about 85 percent) 
occurs at the leading edge of the vehicle. The total ablated weight for 
a 10,000-pound-gross-weight vehicle decelerating from an initial velocity 
of twice circular velocity is estimated to be 5 percent of gross weight. 
Modifying the constant G load flight path by a constant-angle-of-attack 
segment through a flight- to circular-velocity ratio of 1.0 gives essen- 
tially a "point landing" capability but also results in an increased 
total heat input to the vehicle. 


INTRODUCTION 

Studies of interplanetary missions indicate that atmospheric braking 
can be markedly more efficient (in terms of weight requirements) than 



2 


propulsive braking (ref. l) . To apply atmospheric braking at the entry 
velocities associated with interplanetary trips , several problems must 
be studied. One of these is the question of whether the achievable aero- 
dynamic entry corridors (e.g., refs. 2 and 3) are compatible with the 
corridor capability of the planet approach guidance system. A second 
problem is protecting the vehicle from the heat generated in decelerating 
from the entry velocity. 

In terms of acceptable G loads , entry to Earth at velocities 
greater than circular appears feasible. The advent of ablation heat- 
protection systems, in contrast to radiation- and heat-sink systems, also 
makes possible the consideration of supercircular entry speeds from the 
thermodynamic point of view. To minimize the weight of the heat- 
protection system requires the use of "deceleration flight paths" that 
have inherently low total heat inputs. The deceleration flight path is 
defined as that segment of the total flight path in the atmosphere which 
is flown for the purpose of decelerating the vehicle from an initial 
velocity and altitude to a desired terminal condition. It is subsequent 
to the maneuver which determines the entry corridor. The velocity change 
during the entry maneuver may be small (ref. 3), so that to a first ap- 
proximation the initial velocity for the deceleration flight is equal to 
the entry velocity. 

A number of studies of the heat input and flight mechanics of de- 
celeration flight paths have been made (refs. 4 to 7). The stagnation- 
point heat input per unit area nondimensionalized by functions of nose 
radius Rp, lift-drag ratio L/D, wing loading w/A, and drag coefficient 

Cp is the parameter usually presented. Items not considered in these 

studies which are dealt with in the present analyses are: (l) the inter- 

relation among the parameters Rp, L/D, W/A, and Cp required by the 

vehicle geometry, (2) the heat input to the windward side of the vehicle, 
(3) the integration of the heat input per unit area over the vehicle 
area, and (4) the effect of differentiating between the vehicle maximum 
and operating lift-drag ratios. The parameter heat input per unit ve- 
hicle gross weight is used as a measure of the heat input to the vehicle. 
This parameter is closely related to the ablated weight as a fraction of 
the vehicle gross weight. 

The present analysis yields closed-form results. The flight mechan- 
ical approximations made are less restrictive than the assumptions used 
to arrive at existing closed-form solutions. The characteristics of 
the following types of flight paths are derived: (l) constant angle of 

attack, (2) constant G load, (3) constant heat rate, (4) constant alti- 
tude, (5) constant free -stream Reynolds number, and (6) "modulated roll." 
The following constraints are imposed: that the vehicle (l) remain in 

or near the atmosphere, (2) stay within its aerodynamic capabilities, 
and (3) stay within the G load capabilities of the payload and 
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structure. A Reynolds number constraint is also discussed. The total 
heat input to the vehicle and the G load history are discussed as a 
function of the type of flight path, initial velocity, and the initial 
G load. Finally, a method of controlling landing location is discussed. 


ANALYSIS 

The general objective of this analysis is to develop closed-form 
solutions for the characteristics of the deceleration flight path. 


Equations of Motion 

The equations of motion including gravitation and aerodynamic forces 
in the directions tangent and normal to the flight path may he written in 
the nomenclature of figure 1 as 


1 dV _ -CpP^Ar 
g d-r ~ 2m 


sin cp 


-C D pV^Agr 

2W 


- sin cp 


( 1 ) 


V d£ C LP v2Ar ,r=2 

g dr 2m ' 


- l)cos cp = 


C L pV^ Agr 
2W 


+ (V‘ 


- l)cos cp 


( 2 ) 


The symbols are defined in appendix A. 

The following assumptions are made: (l) that the flight -path angle 

relative to the local horizontal is always small, so that cos cp is ap- 
proximately unity and sin cp is approximately cp, and ( 2 ) that sin cp 
is small compared with the nondimensional term representing the aerody- 
namic forces in equation (l) . 


Equations (l) and (2) then reduce to 


1 dV = _ Q S 
g dx 


V d£ 

g dt 



1 


(3) 

(4a) 


where G is defined as 


G = 


CppV^Ar 


CRpV^Agr 


2m 


2W 


(5) 
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G represents the net force on the vehicle in terms of the planet g f s, 
or is the net acceleration vector felt by the pilot in planet g T s. 

For future reference, equation (4a), with the aid of equation (5) 
and the definition of (appendix A), may also be written as 


_L _ 

mg 



+ I *2. 

g dr 


pV 2 = 


2m (l ~ V 2 + ~ 2w(l - V 2 + ^ 


C L Ar 


C L A gr 


( 4 b) 


(4c) 


Also from figure 1, the component of lift in the vertical direction is 
related to the total aerodynamic lift a and the roll angle 0 by 

C L = C L,a cos 6 ( 6 ) 


Vehicle Aerodynamics 

The vehicle aerodynamics (fig. 2) are calculated using the equations 
and assumptions given in appendix B. The vehicle design is specified by, 
or characterized by, its maximum lift-drag ratio (L/D) max - In general, 

for a given lift coefficient these curves have double values of C-q, C^, 
and a. In the present study the low range of angle of attack is used. 


Flight Paths 

Many possible flight paths can be considered; the present report 
considers six particular flight paths and a generalized flight path. The 
flight paths are defined as follows: 


(1) 

fig. 2) . 

Constant 

a: 

0=0; , C]}/C^, and Cp/Cp are constant 

(2) 

Constant 

h: 

0 = 0; p is constant. 

(3) 

Constant 

G: 

0 = 0; CrpV 2 is constant (see eq. (5)). 

(4) 

Constant 

Re : 

0 = 0; pV is constant (see eq. (17)). 

(5) 

Constant 

4 = 

l/2— 3 

0 = 0; p ' V is constant (see eq. (10)). 
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(6) Modulated roll: Q is varied from 0° to 90° (or 180°) in a man- 

ner to satisfy equations (4) and (6); the vehicle angle of attack is de- 
fined to he constant, and hence > c d/ c R> and C L,a/ C R are constant. 


Heat Transfer 

The heat transfer for the present analysis is estimated assuming a 
laminar boundary layer. References 8 and 9 present a comprehensive re- 
view of heat-transfer theories. These have been found adequate up to 
Earth circular speeds and are assumed herein to be applicable also at 
supercircular speeds . 

Reference 10 gives the following form for stagnation-point heat 
transfer : 


15.5X10 


" 9 X2j/ 2 p 1 / 2 V 3 

"W 5 



(7a) 


where j = 0 for the two-dimensional leading edge considered in the 
present analysis (j = 1 for a three-dimensional or spherical leading 


edge) . For entry from high velocities 
(7a) in a form using V, 



« 1.0 


and with equation 


a = 15.5X10- 9 (gr) 3 / 2 


R 


1/2 


■N 


(71) 


This equation may be modified to include the effect of sweep and nose 
radius angle by the relation suggested in reference 8 (pp. 19 and 22): 


— = cos v A c cos A = JjL - sin 2 A cos 2 aj V/ ^ cos A 


Combining equations (7) and (8) gives 


( 8 ) 




= 15.5X10 


-9 


(gr) 3 / 2 


3 l/2y 3 

w 


cos 


V A„ 


cos 


A 


(9) 


The theory and experiments of reference 11 indicate that, for the high 
velocities of the present study, v has a value of 1.5 for sweep angles 
from 0° to 60° and becomes less than 1.5 for sweeps greater than 60°. 
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(The theory and experiment of ref. 12 indicate that the value of v may 
he as low as unity.) 

The symbols defining the vehicle geometry and various flow compo- 
nents are illustrated in figure 3. In calculating the heating rates and 
total heat inputs, a family of vehicles was considered which are delta- 
wing configurations having a wing sweepback angle A, a constant thick- 
ness d, and a constant nose radius Rjj. The configurations have no 

separate fuselage, and the volume in the wing is considered the usable 
volume . 


The following analogously simple equation for the heat transfer to 
the windward side of the vehicle is developed in appendix C: 


4 W = 15 .0X1C> k n gr 


JHE 

T Ax I c 


cos 9 1 


= 10.2 k. 




G f Ax I cos 01 


v 2 


( 10 ) 


This result is independent of angle of attack between about 3° and 30° 
for a constant roll angle 6. 


General Form of Equations for Flight-Path Characteristics 

In addition to those characteristics involved directly with the 
equations of motion there are other flight-path characteristics of in- 
terest and these are given in general form here. 

Altitude . - With the assumption of an exponential atmosphere, the 
flight-path altitude may be calculated if the flight-path density is 
known : 


(11) 

For Earth the values 3" 1 = 2.35X10 4 feet and p gL = 0.0027 (slug) (ft -3 ) 

are used as suggested in reference 4. Also for later use it is noted 
that 


dh _ JL _ 
dp Pp 


( 12 ) 


Radiation equilibrium temperature . - With the flight path and ve- 
hicle design known, equations (9) and (10) give the heating rate per 
square foot q. In some applications the radiation equilibrium tempera- 
ture, determined by the heat rate, is of interest; it is given by 


^ = 



1/4 


( 13 ) 
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Free-stream Reynolds number . - With the flight -path density p 
known as a function of V and the vehicle characteristic length, the 
free-stream Reynolds number may be calculated from 

Re = ^ VI? (14) 

p 


By use of equation (5), equation (14) may alternatively be written as 


Re , 

pVs r VC R 


(15) 


The free-stream viscosity is assumed constant and independent of altitude 
for the present analysis, p = 0.37X10® slugs per foot-second, and x is 
taken as f plus the root chord c (see fig. 3). 

Local Reynolds number . - Local Reynolds number is related to the 
free-stream Reynolds number by a function of the free-stream Mach number, 
the wing-sweep angle, and the angle of attack. The method of calculation 
of the ratio of local to free-stream Reynolds number Re^/Re is given 

in appendix D and assumes y = 1.2. The local Reynolds number is then 



The value of Re^/Re depends on the value selected for y, and more 

generally on the real gas properties. A more precise determination of 
the local Reynolds number than that of appendix D is of little value at 
present, because the transition Reynolds number is not accurately known. 

Flight-path angle . - As mentioned previously, it is assumed that 
the flight-path angles are small, and hence 

V 

sin cp ~ cp ~ (17) 


But 


V dh dh dp d V 
v " dx “ dp dV dr 


and by use of equations (3), (5), and (12) 


C D S V dp 

9 = 77 — 

2(3 J 


(ie) 


( 19 ) 


The preceding characteristics are point functions along the flight 
path. The succeeding characteristics result from an integration along 
the flight path. 

Time of flight . - The time of flight is obtained by integrating the 
time increment from equation (3) as follows: 



Range . - Using equation (3) for dx and approximating cos cp by 
unity give the range traversed during the deceleration: 



Total heat input . - One of the important criteria of merit for a 
heat -protection system is its weight as a fraction of the vehicle gross 
weight. If, for example, the heat-protection system is of the ablation 
type and has a characteristic heat-absorbing capacity C in Btu per 
pound, then the ablated weight may be expressed as 


i i a 

W C ¥ 


( 22 ) 


The present analysis deals with the parameter Q/W which results from 
a double integration (along the flight path and over the vehicle). The 
integration along the flight path yields the total heat input per square 
foot, which is given by 


1 = 






(23) 
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where q is given by equations (9) and (10) and dr is again given by 
equation (3) . 

The integration over the vehicle has the general form 

M ff di i + 1 ff ^ dA dZ (24) 

*'A W A z 

This integration is discussed in detail in appendix E. 


Sample Development of Flight-Path Characteristics 
for a Constant Reynolds Number Path 

The case of "constant free-stream Reynolds number" is considered 
here in detail. The velocity ratio V is taken as the independent vari- 
able. The free-stream Reynolds number is given by equation (14) or (15) , 
if the vehicle is specified (i.£., f+c and w/A are constant and given) 
Then for Re to be constant, pV must be constant. Since G (eq. (5)) 
and all the terms collected on the left side of the equation 

^ - fe} 


are defined as constants, the right side must also be a constant. This 
is the relation given in column 2 of table I. The first column identifie 
the type of flight path. The term in braces is a constant for the flight 
path and is specified if the path Reynolds number is specified. 


Equation (25) may then be solved for p as a function of V to 
give column 3 in table I: 


P = 


2 £ 1 f G- 1 1 
A gr 


(26) 


Because {g/C^v} is a constant along the flight path, 



* 


( 27 ) 
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which is the result given in column 4. To determine the flight-path 
angle cp, the derivative dp/dV is required in equation (19) . This may- 
be obtained from equation (26) : 


d£ = _ W 1 f g1 _1 
dV ” A gr |C R V[V2 


(28) 


Combining equations (19) and (28) leads to the result in column 5: 



cp = - 


With cp having been determined, the term ^ ^ can be discussed, 
equation (29) 


V d£ C D f G 

i iav 

g dr ~ 3r 

v e dT 


(29) 


From 


(30) 


where the term — ~ is given by equation (3), and G in that equation 

is given by equation (27) (or column 4). With these substitutions equa- 
tion (30) becomes 


V deg _ _1_/g1 2 /9d' \ 2 ^2 = 

g dr " (3r jO^Vj \C R j 


(31) 


The term ^ ^ (column 6) is equal to -k^ e (column 19) . For 
and Cft constant, an assumption required later, k^ e is constant. 

Equation (4a) may now be solved for the lift coefficient, using the 
result of equation (3l) and equation (27) for G, to arrive at the re- 
sult in column 7 : 



The free-stream Reynolds number (column 8) is given by equation (14) 
and is a constant for thi_s case. More generally it may be obtained from 
the relations of p to V like equation (26) or column 3 of table I. 
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Obtaining the local Reynolds number (column 9) requires knowledge 
of the wing sweep angle and the angle of attack. The wing sweep is part 
of the vehicle design specification. The angle of attack is determined 
by figure 2 using Cp from equation (32) and the vehicle (L/D) max . 

The heating rate to the leading edge q^ (eq. (9)) now depends only 

on the parameter p, which has been determined, and on the independent 
variable V. This result is shown in column 10. The heating rate to 
the windward side q w (eq. (10)) depends on L/A, which may be determined 

from equation (4b) and the result of equation (3l) or column 6 as fol- 
lows: 


| = | (1 - V 2 - kRe) (33) 

Making this substitution in equation (10) gives the result listed in 
column 11 of table I. 

The point functions have now all been determined, and it remains 
to determine the integrated functions. The time of flight may be ob- 
tained from equation (20) and the equation for G from column 4 (or 
eq. (27)), where the term in braces is a constant, so that 



The variation of and (or just Cp) as a function of V re- 

quired in equation (34) can be obtained from the curves of figure 2 and 
equation (32) for Cp. Hence, equation (34) can be integrated graphi- 
cally or by numerical methods for the most general cases. The present 
report, however, deals only with the more simple cases where Cp/Cp and 

Cp can be approximated by constant values and closed-form solutions can 

be thus obtained. This assumption is different from the assumptions for 
the case of constant a for which Cp/Cp in addition to Cp/Cp and 

Cp is constant. With the assumptions of Cp/Cp and c R constant. 
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which are noted in column 12, equation (34) becomes 



which integrates to the result in column 13. 

By similar substitutions and assumptions, equation (21) for this 
range becomes 



which integrates to the result in column 14. 


Again with the similar substitutions and approximations, and by 
use of the result in column 10 or equation (9) for q^, equation (23) 

for q^ becomes 



which integrates to the result in column 15. 

When column 11 or equation (10) with 0=0 is used for q^ in 
equation (23) , the integral equation becomes 


q w =!0.2 k c V“ 




_ 0 _ 

C R V| 





V 2 - 


s Re I 


V dV 


(38) 


which integrates to the result in column 16. 

The integration over the vehicle to obtain Q/w indicated by equa- 
tion (24) is carried out in detail in appendix E and gives the results 
shown in columns 17 and 18. Alternative forms for the denominator of 
columns 17 to 18 are shown in the bottom row. 
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The results for all the other flight paths except constant a can 
he derived in a similar manner. The case of constant angle of attack is 
discussed in appendix F. 


RESULTS AND DISCUSSION 

The results of the analysis for all the flight paths studied are 
summarized in tables 1 and II. 


Tables of Equations 

The first column of table I gives the name of the flight path. The 
term in braces in the second and subsequent columns is constant over the 
flight path. The equations with the exception of those in columns 10, 

11, and 15 to 18 apply to any planet. The results in the aforementioned 
columns are given in table II in a form which is applicable to any planet 
that has an atmosphere with heat-transfer characteristics simila£ to 
Earth. The equations for G- are in terms of planet g T s, and V is the 
flight velocity divided by the planet circular velocity. 

Three heating parameters are presented: 

(1) Heating rate q (inBtu/(sq ft)(sec)). This is the parameter 
of interest, for example, for radiation cooling. 

(2) Heat input per square foot of vehicle q, which depends on the 
location with respect to the stagnation point of the area being con- 
sidered. (The stagnation point is not a fixed geometric point except 
for the cases of constant a and modulated roll.) The parameter q re- 
sults from the integration over the flight path of the heat input to a 
given location. 

(3) Heat input Q = (Q^ + Q^), which is related to the ablated weight 

and involves an integration over the vehicle surface as well as over the 
flight path. 

The heating parameters are presented for two regions on the airplane: 
(l) the region in the vicinity of the leading edge or stagnation point, 
where the heating is a function of nose radius R^, nose radius coordi- 
nate angle A, and sweep angle A, and (2) the windward side of the ve- 
hicle, where the heat input is a function of the distance x from the 
stagnation point. The heat input to base and lee regions is not con- 
sidered. 
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The term (column 17) accounts for the reduction with angle of 

attack of both the effective sweepback angle and the area of the leading 
edge subject to heating. The value of k^ is approximately unity. The 

term k c (columns 11 , 16, and 18) accounts for the relative lift generated 
by the windward side and the leading edge. The term k w accounts for 
the length of boundary layer ahead of the windward side. The term k 
with a subscript to designate the flight path (e.g., % e ) accounts for 

the amount that the flight -path curvature differs from the planet surface 
curvature. Both k c and approach unity for slender (or high 

(L/D) max ) configurations and have a value of about 0.5 for the 

(L/D) max = 1.0 configuration of this report. 

In table I the heading "Modulated angle of attack" indicates the 
flight paths for which the lift vector is constrained to the vertical 
direction. For all but the constant a case, the assumption of con- 
stant C;d/Cr restricts the angle of attack to the vicinity of 0° or 90°, 

as may be seen from figure 2. The present results are for a near 0°. 

"Modulated roll” in table I designates a specific case where angle 
of attack with respect to the free-stream direction is constant (hence, 

°R> c d/°R^ and C L,a/ C R are constant) and the variation of the roll angle 

0 with _V is given by the equation in column 1^ This relation between 

0 and V gives a constant G and 0 = 0 at Vj_. 

Comparison with Other Analyses 

To obtain closed-form solutions, a number of simplifying assumptions 
have been made. As an indication of the accuracy of the present analysis 
in the stagnation region, a comparison is made for the constant-angle- 
of-attack path with the results from machine integrations of the equa- 
tions of motion and heating presented in reference 4. The procedure for 
reducing the results of the present report to the form of reference 4 is 
given in appendix G. Figure 4 presents the horizontal component of the 
G load Gjf and the heating-rate parameter q as a function of velocity 

ratio V for an L/D of 0.1 and 1.0 and an initial velocity ratio V 
of 0^99. The heat-input parameter q for decelerating from V = 0^99 
to V = 0 is given as a function of L/D. In the calculation of c[, 
the function F of table I was evaluated at the maximum value of V. 

The present analysis generally gives the greater heat input. In all 
cases the agreement between the two analyses is considered satisfactory. 
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To indicate the accuracy of the present analysis in estimating the 
heat input to the windward side of the vehicle, a comparison is made with 
the results of reference 13 for the constant-angle-of -attack path. The 
results of reference 13 have been extrapolated to the condition of zero 
wall temperature and apply for the case of F = 1 in table I . Figure 5 
compares the heating-rate parameter q _(x/(w/a))V^ and the heat input 


per unit span parameter 


/ 


q dx 


/c'W 

\ A . 


1/2 


as a function of velocity ratio V 


) 

for an initial velocity V of 1.0. Also shown is the variation of heat- 
ing rate with angle of attack for constant lift. In general, the agree- 
ment between the two analyses is satisfactory. Compared with reference 
13, the present analysis indicates increasingly higher heat inputs with 
increasing angle of attack. The maximum angle of attack to which the 
present analysis is applicable depends on the degree of approximation 
considered acceptable. The present analysis is about 15 percent higher 
at a = 30°. 


In general, the accuracy of the present analysis is adequate for 
preliminary studies. 


Some General Considerations 

For atmospheric braking it is desired to determine a combination 
of vehicle geometry and flight path which will result in a low heat in- 
put to the vehicle. Generally there is a substantial increase in con- 
vective heat input to the vehicle if the boundary layer on the vehicle 
is turbulent rather than laminar. Thus it is desirable to restrict the 
flight path to conditions for which laminar boundary layer can be main- 
tained. An indication of the possible nature of the boundary layer is 
given by the magnitude of the Reynolds number. Two Reynolds numbers can 
be considered, the free-stream value, and a local value at the edge of 
the boundary layer. The free-stream Reynolds number depends on the 
vehicle design parameters w/A and c and on the flight-path parameter 
G (see eq. (17) or (18)). The local Reynolds number is further influ- 
enced by the vehicle design factors of leading-edge bluntness and sweep 
and the operating parameter of angle of attack. Both local and free- 
stream Reynolds number are presented in the section Point Functions. 

The local Reynolds number is interrelated with the problem of hot- 
gas radiation. At entry velocity ratios V greater than 1.0, the heat 
input from hot-gas radiation can equal or become larger than the con- 
vective heat input to an unswept stagnation region. The lowest heat in- 
put to the vehicle will usually exist if hot-gas radiation is avoided 
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where possible. This requires that the vehicle surface angles be small 
with respect to the free-stream direction to avoid high air static tem- 
peratures. In general, unsweeping the wing leading edge decreases the 
local Reynolds number but tends to aggravate the problem of hot-gas radia- 
tion. The choice of wing sweep angle is, therefore, a compromise. The 
calculations for the present study were made for a sweep angle of 60°. 

Also to avoid hot-gas radiation the low-angle-of-attaek rather than high- 
angle -of -attack region of aerodynamic operation was chosen. The choice 
of the low-angle -of -at tack range and the 60° sweep is assumed to reduce 
the heat input from hot-gas radiation to a low value compared with the 
convective heat input. Only the convective heat input is considered 
subsequently in this report. 


Point Functions 

The evaluation of the various flight paths for decelerating in the 
Earth atmosphere is discussed with the aid of figures 6 to 11 in which 
is plotted the variation of the point functions of table I with V. A 
second abscissa shows elapsed time in minutes with time equal to zero at 
V = 2.0. To arrive at numerical results, a vehicle design defined by 
= 20 pounds per cubic foot, W = 10,000 pounds, A = 60°, and 

(L/D) max = 1.0 was assumed. (Higher values of (L/D) max result in 

larger heat inputs.) For a vehicle like that of figure 3 and the above 
numerical values, the following vehicle characteristics may be determined 
by the equations given in appendix E: w/a = 83 pounds per square foot, 

R^ = 2.06 feet, and c = 14.5 feet. A length of 18.8 feet was used in 

the calculation of the Reynolds numbers. Vehicle operation is in the 
low-angle-of-attack range, and the assumptions are made that C-p = -Cp Q 

and C-q/C^ = 1.0, although for some cases more precise approaches to the 

calculation of the point functions are possible. All the flight paths 
have a common point, at V = 2.0, G = 10 g*s, and 0 = 0. This results 
in approximately the same values of a and h, about 5.3 and 162,000 
feet, respectively. When appropriate, the following physical limits are 
noted on the curves: |Cp^ ma ^;] = 0.9, G max = 10 g's, b min =0, and 

q = 190 Btu/(ft~ 2 ) (sec) , or t— = 4460° R for an emissivity of unity. 

Flight -path angle . - The fact that for all cases the flight-path 
angle cp is small except in the vicinity of V = 0 justifies the assump- 
tion of small angles. For all cases except constant a (fig. 6), cp is 
also continuous. For the_constant a case, cp changes dis continuously 
from plus to minus near V = 1.0 and the altitude becomes indefinite. 

This situation is physically approached by a shallow skip out of the 
atmosphere at V = 1.0. In spite of this indefiniteness in the altitude, 
the heat input to the vehicle remains finite as is shown later (see also 
columns 17 and 18 of table i) . 
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Altitude and lift coefficients . - For all cas£s except constant h 
(fig. 7) the altitude becomes less”” than zero for V approaching zero. 

On the constant h path excessive lift-coefficient requirements are en- 
countered at V l s approaching zero. The constant^ Re path (fig. 8) 
also requires excessive lift coefficients at low V l s. In a practical 
flight these and other low-speed anomalies are usually precluded by 
specifying some terminal condition of velocity, lift coefficient, or 
altitude . The actual flight path may deviate from the specified path in 
the low- speed range with little effect on the total heat input calculated 
from an initially supercircular velocity, as is shown later. 

The lift coefficient is negative; that is, lift force is directed 
toward the Earth for flight velocities greater than circular. This 
occurs fundamentally because the centrifugal force tending to "throw" 
the vehicle out of the atmosphere is greater than the gravity force and 
the difference in these forces must be counteracted by the lift. For 
all cases except constant G (fig. 9), constant q (fig* 10) , and modu- 
lated roll (fig. 11), the lift coefficient is positive (directed away 
from the Earth) for subcircular velocities, because at subcircular speed 
the centrifugal force is less than the gravity force. However, to satisfy 
the assigned values of q and G for the cases of constant q, constant 
G, and modulated roll, negative lift coefficients are required at sub- 
circular speeds as well as at supercircular speeds (because of the term 
(v/g)(dcp/dt) in eq. ( 4 )). 

G load . - For the constant G path (G = 10 g's), the abscissa of 
figure 9 shows that the time of deceleration in the region of strong 
heating, that is, from V = 2.0 to V~0.5, is about 2 minutes. This 
combination of g's and time is about the presently accepted limit of 
human tolerance if useful functions are to be performed by the man. For 
this reason G max = 10 g's was used in the present examples. 

For all cases except constant < 3 . (fig. 10) the G load is at or 
below the initial G value at V =_2.0. For the case of constant q 
the G increases with decreasing V to maintain q constant. At a 
V higher than 2 the becomes excessive. These same characteristics 

exist if other values of constant 4 are specified. Thus it is con- 
cluded that the constant q path does not have application over a wide 
range of V. 

Angle of attack . - The angle of attack is shown as positive. For 
the flight paths where the lift coefficient varies from minus to plus 
along the flight path, one has the choice of permitting the sign of the 
angle of attack, to change or of rolling the vehicle 180°. In the latter 
instance the same side of the vehicle remains the windward side. The 
vehicle angle of attack is for all cases equal to or less than the angle 
of attack at V = 2.0, where a is about 5.3° (except for the very low 
V's). From figure 2 it can be seen that Oq, C^/C^, and vary less 


than ±5 percent from a mean value over the angle of attack range 0° to 
5.3°. More generally, for = 10 g's and ( Cp/ t = (l - V 2 )/G i 

^eq. (4a) with ^ = 0 j , C D , for example, varies less than ±7 percent 
from a mean value up to = 2.25 for (L/D) max = 1 and 2, the values 
of (L/D) max examined. An increase in G^ will increase Vj_, for which 

the variation of the coefficients is small. Hence, for many problems of 
interest it is a good approximation to use constant values for Gq, 

Gp/C^, and C R as is required to arrive at some of the integrated func- 
tions . 

Free-stream Reynolds number . - The free -stream Reynolds number Re 
is also shown in figures 6 to 11. The case of constant h shows a 
linearly decreasing Reynolds number with decreasing V. The case of con- 
stant G for G = 10 g's shows increasing Reynolds number with decreasing 

V and values much larger than for the other cases in the lower V range. 
The various flight paths illustrated thus offer a choice of increasing, 
constant, or decreasing Reynolds number with decreasing V. 

Local Reynolds number . - The transition from laminar to turbulent 
flow in the boundary layer is more significantly related to the local 
Reynolds number than to the free-stream Reynolds number. The variation 
of the local Reynolds number Re^ with V is shown for the condition 

that the wing leading edge is blunt and swept 60°. All the flight paths 
except modulated roll (fig. 11) show a dip in Re^ in the vicinity of 

V = 1 followed by a rise in Re^ for V <1.0. For the flight paths of 

constant G and constant Re,_this dip is associated with a lower angle 
of attack in the vicinity of V = 1.0. For the constant a case it is 
associated with a low air density (see fig. 6). In the vicinity of 

V = 1.0, the local Reynolds number for both the constant a (a ~ 5.3°) 
and constant G (G = 10) cases is substantially lower than for the 
modulated roll case (G = 10 g's and a = 5.3°). This illustrates a dis- 
advantage of using modulated roll. The high local Reynolds number for 
modulated roll is associated with the combination of the high air den- 
sity and the 5.3° angle of attack at V = 1.0. The effect of angle of 
attack on the ratio Re^/Re is discussed in appendix D. In general, 

for the flight paths considered, it appears less difficult to achieve 
low local Reynolds number at supercircular than at subcircular speeds. 

The constant a, constant h, constant Re, and constant G flight 
paths show rather small variation in Re^ with V and little variation 

among the several flight paths over the speed range from V = 3.0 to 
0.8. In this speed range the general level of the local Reynolds number, 
0.7X10 6 , suggests an extensive run of laminar boundary layer on the ve- 
hicle is possible. 
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Heating rate . - The values of the heating rate q shown are for a 
vehicle density = 20 (lb)(ft“3) and apply to the swept leading edge 

(A = 60°) (note that, in column 10 of table I, W/R-^A = 2p^) . For all 

cases except constant q, values of q increase with increasing V. 

The temperature limits of presently known materials (e.g., ceramics) are 
about 4460° R. By equation (16) this corresponds to 

q = 190 (Btu) (ft - ^) (sec - -*-) . For this value the best case, constant a, 
shows that entry velocity ratios V up to 1.3 may be possible with 
radiation cooling, for the present vehicle and flight path. (The flight 
path and vehicle design shown do not yield the lowest possible radiation 
q temperatures.) For entry velocity ratios greater than approximately 1.4 

radiation cooling is probably not possible, and hence there is interest 
in the total heat input to the vehicle . 

Maximum entry velocity . - From equation 4(a) for = 10 g's and 
the maximum values of C^/C^ shown in figure 2, it can be determined 

that entry velocities of 2.8 times circular are aerodynamically possible 
with a vehicle with an (L/D) max Still higher entry velocities 

are possible with a vehicle with a higher (L/D) max . 

Summary of point functions . - The more important conclusions from 
the discussion of the point functions may be summarized as follows: 

(1) All the flight paths considered except constant q are or ap- 
proximate physically possible flight paths. 

(2) For G ± = 10 g's and (L/D) max = 1.0 or 2.0 the assumption of 

Cr, Ojj, and. constant is a valid approximation up to an entry 

velocity ratio V of about 2.25. 

(3) With the vehicle designs and flight paths assumed, radiation 
cooling appears impossible for entries at supercircular speeds corre- 
sponding to V > 1.4. Hence, the total heat input into the vehicle is 
important . 

(4) The local Reynolds number for many of the flight paths of in- 
terest is about 0.7X10 6 . This value suggests the possibility of a 
laminar boundary layer. 


Total Heat Input 

For the assumptions of laminar boundary layer and hot-gas radiation 
negligible because of low vehicle surface angles and with heat input to 
base and lee regions neglected, the total heat input to the vehicle is 



given "by the sum of columns 17 and 18 of table I or II . The total heat 
input depends on many factors, and these may be categorized as follows: 

(1) Factors that depend on the flight plan:_ terms in braces or in 
column 2 that define the variation of G with V and the terms in 
brackets in columns 15 and 16 that depend on the initial and final veloc- 
ity of the deceleration flight path. 

(2) Factors that depend on the vehicle design: vehicle density, 

p f ; gross weight, W; wing loading, W/A; and wing sweep angle, A. 

( 3 ) Factors that depend on both vehicle design and flight plan: C ^ 

and L/D or C^/C^ and C^. These depend on design through the con- 
figuration maximum lift -drag ratio and on the flight plan through the 
operating lift-drag ratio or angle of attack. 

The heat input along families of flight plans related in two ways 
to the curves of figures 6 to 11 will be considered. For the first family 
the flight plans start at various points along the curves of figures 6 to 
11 and the flight path then follows these curves. The second family is 
defined by the condition that, at the initial velocity V^, the initial 

G load G^ is 10 g's. 

Figure 12 considers the first family of flight paths, which have in 
common a net acceleration of 10 g's at V = 2.0, and presents the heat 
to the windward side of the vehicle and the total heat input Q 

(the sum of the heat inputs to the windward side and to the leading-edge 
region) as a function of Vj_ . Though not shown by the curves, from the 
equation of tables I and II, the heat input Q equals zero for Vq = 0. 

Each point in figure 12 can be interpreted as giving the heat input 
for a complete "deceleration flight path" having a terminal velocity of 
zero and the initial velocity and initial G given by figures 6 to 11. 

The history of the flight is also given by figures 6 to 11. For the 
supercircular speed rangp the flight paths in order of increasing total 
heat input for a given Vj_ rank as follows: constant G, modulated 

roll, constant Re, constant h, and constant a. At = 2, the con- 
stant a case has about twice the heat input of the constant G case. 

The curves of figure 12 can also be interpreted to show the distri- 
bution of heat input as a function of V. For example, for a flight 
starting at any initial velocity, the contribution to the heat input in 
decelerating from V = 0.5 to 0 is given by the Q at V = 0.5. With 
this interpretation for all the flight plans the heat input for decelerat- 
ing from V = 0.5 to 0 is generally an order of magnitude less than the 
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total heat input for decelerating from any super circular speed. A precise 
determination of the total heat input at the low speeds is not essential 
to an accurate estimate of the total heat input in decelerating from 
supercircular speeds. This implies that neglecting the low-speed anom- 
alies previously discussed in regard to the point functions has little 
effect on the accuracy of calculating the total heat input for the de- 
celeration from supercircular speeds, and therefore these anomalies have 
been neglected in the present calculations. For most flight paths the 
total heat input to the windward side of the vehicle is about 15 percent 
of the total heat input, for a vehicle with an (L/D) max of 1.0. 

The second family of deceleration flight paths considered is shown 
in figure 13. Each point in the figure represents a complete decelera- 
tion flight path having the given by the abscissa and = 10 g*s. 

The values for = 2.0 for all the flight plans and for the complete 

curve for constant G are the same for figures 12 and 13. In the super- 
circular speed range the flight paths generally rank in the same order 
as for the preceding family of curves. The equations for the constant a 
path are of questionable accuracy for the very low L/D's (l/D < 0.03) 
that result for close to unity, and hence results are not shown at 

V less than 1.2. 

The relation of the ablated weight to the total heat input is given 
by equation (26). An ablation system based on the extrapolation of the 
data of references 14 and 15 to the conditions at supercircular speeds 
and a system utilizing the boiling of lithium (ref. 16) both indicate 
the possibility of an effective heat-absorbing capacity of about 10,000 
Btu per pound of ablated weight. Using this value for C in equation 
(26) for = 2.0 and G^ = 10 g*s gives values of the ablated weight 

as a fraction of the vehicle gross weight of 0.039 and 0.085 for the 
flight paths of constant G and constant a, respectively. Halving the 
value of C would double the ablated weight . From the consideration of 
the ablated weight it is concluded that entries at V = 2.0 and greater 
are feasible . 

The distribution of heat between the windward side and the leading 
edge of the vehicle depends on the G load, as illustrated in figure 14 
for a family of constant a flight paths with = 2.0 and for a ve- 
hicle with an (L/D) max of 1.0. As the initial G load decreases, the 

vehicle angle of attack and initial flight altitude must increase to 
satisfy equations (4a) and (5). Because of the assumption of constant 
and Cp/Cp the other flight plans cannot be considered over this 

wide a range of G or a. The curves show that the heat input to both 
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the leading-edge region and the windward side of the vehicle increases 
with decreasing G load. Reducing the initial G^ from 10 to 4.25 g's 

about doubles the total heat input. The value G i = 4.25 g’s is the 

minimum at which the constant a path can be initiated. Also the pro- 
portion of the heat input to the windward side of the vehicle increases 
from about 19 percent at = 10 g's to 47 percent at Gj_ = 4 . 25 g 1 s . 

The increased angle of attack accounts for the increased proportion of 
heat to the windward side of the vehicle. 

While the heat input for the constant a, G^_ = 4.25 g path is 2.2 
times that for the constant a* G^ = 10 g path, it is 4.9 times larger 

than that for the constant G = 10 g path. This indicates the importance 
of both the G level and the type of deceleration flight path to the 
vehicle heat input. 

Figure 14 also shows the Reynolds numbers at = 2.0. Decreasing 

the G load decreases the free-stream Reynolds number almost directly, 
as would be expected from equation (18) (for constant C R ) . Increasing 

angle of attack, however, increases the ratio of local to free-stream 
Reynolds number (see appendix D), so that decreasing Gj_ has essentially 

no effect on the local Reynolds number. 


Control of Landing Location 

One of the practical requirements that may be placed on an entry 
vehicle is that it have the capability of landing at a predesignated 
location. The problem of achieving a low total heat input on a flight 
path that will also yield control of the landing location can be studied 
to a first-order approximation with the present solutions. 

To land at a predesignated location requires control of motion both 
down and cross range with respect to the Earth 1 s surface. By column 14 
of table I the down range for the constant a case is indefinitely large. 
For this case the analysis implies that at V = 1.0 the vehicle glides 
up out of the atmosphere into an orbit of indefinite duration. (The de- 
tails of this maneuver are not defined by the equation. If the vehicle 
exits from the atmosphere at a small angle as indicated by fig. 6, a 
small rocket may be required to "circularize" the orbit, and then a second 
small rocket may be required to initiate reentry.) Cross range movement 
may be obtained by waiting for the Earth to rotate through the plane of 
the orbit If, for example, the vehicle had glided into a polar circular 
orbit, all of Earth's surface passes through the plane of the vehicle in 
l/2 revolution of Earth or in 12 hours. This time corresponds to about 
8 revolutions of the vehicle about Earth. A down-range capability of 8 
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revolutions thus yields a control of cross range. With the aerodynamic 
cross-range capability such that half the distance between the traces on 
the Earth l s surface of two successive orbits can be traversed during the 
deceleration from the circular orbit , the vehicle can potentially land 
at any predesignated location. 

The present discussion deads with the down range and total heat in- 
put for some possible flight plans that_include control of landing loca- 
tion by control of the flight through V = 1.0. The flight plans con- 
sidered are illustrated in figure 15 in terms of the G history against 
V. The flight paths are composites of the constant G and constant a 
flight paths, and all begin at V = 2.0 with an initial deceleration of 
G = 10 g's. Figure 15(a) shows the constant a case modified by a con- 
stant G segment through V = 1.0. The magnitude of the deceleration 
during this constant G segment is designated G m j_ n , because it is the 

minimum deceleration that occurs during the deceleration. Then G^pn 
can be considered an independent variable. For G m j_ n = 0 the general 
path reduces to the constant a case and for G^in = 10 g's to the con- 
stant G case . 

The variation of Q and range with G min is presented in figure 
16. The minimum Q and range occur for G m ^ n = 10, and the values are 

390 Btu per pound of gross weight and 0.032 of the Earth circumference, 
respectively. The maximum Q and range occur at G m j_ n =0. At 

G m i n = 0 the value of Q, is a little more than twice the minimum, and 

the range is indefinitely large. These data show that increased range 
results in an increase in total heat input and that the maximum total 
heat input is finite and occurs for an indefinitely large range. 

Figure 15(b) illustrates a family of flight paths, all of which, 
according to the present^analysis, have indefinitely large range. The 
deceleration begins at V = 2.0 along a constant G path. At a speci- 
fied velocity designated V a the flight path is changed to a constant a 

path, which yields an indefinitely large range. With V a considered an 
independent variable, for V a = 2.0 this case reduces to the constant a 

case previously considered. Figure 17 presents the total heat input 
plotted against V a . Decreasing V a reduces the total heat input. At 

V a = 1.3 the total heat input is 32 percent more than the value for a 
constant G of 10 g's, compared with 120 percent more for V a = 2.0, 
which is the constant a case. 

To obtain control of the range the flight paths of figure 15(b), 
like those of figure 15(a), can also be modified by a constant G segment 
through V = 1 as shown in figure 15(c). From the discussion of the 
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flight paths of figure 15(a), it can be concluded that the heat inputs 
for the flight paths of figure 15(c) will fall between the two curves of 
figure 17 . 

Several comments about the limitations of the present range control 
calculations may be made. The composite flight paths considered in this 
section have mathematical discontinuities in flight-path angle at the 
junctions of the various types of flight paths, so that some transition 
maneuver is required. The details of the transitions, however, are not 
expected to have a significant effect on the total heat input to the 
vehicle . 

The details of the flight path through V = 1.0 for the constant 
a flight plan and the determination of a minimum permissible value of 
V a both require a more detailed study than can be made with the present 

assumptions . 


Generalizing the Results with Respect to Vehicle 
Gross Weight and Vehicle Density 

The numerical results for the heat input have been calculated for 
a specific vehicle with a gross weight of 10,000 pounds and a vehicle 
density of 20 pounds per cubic foot. With the assumption of a laminar 
boundary layer for all vehicles, these results can be applied to vehicles 
of different gross weight and average density by the following reasoning. 
Consider the results of columns 17 and 18 with the substitutions given at 
the bottom of the table. In the present application the product 
( c d/ c r) Cr is approximately Cj^ 0 . For configurations like that of fig- 
ure 3, the drag coefficient at a = 0, C D Q , depends only on the vehicle 

proportions, that is, on the thickness ratio d/c and sweep angle A, 
if the contribution of friction is neglected. The contribution of fric- 
tion is small for the low L/D configuration considered in this report. 
The Q/W for a configuration with (L/D) max = 1 but with any given 
weight and fuselage density may hence be estimated from the curves of 
this report by the relation: 




curves 


x 
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CONCLUSIONS 

An approximate analysis of flight paths for the descent from cir- 
cular and supercircular speeds has yielded closed-form solutions for the 
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following flight-path characteristics: altitude, lift coefficient, G 

load, path angle, free-stream and local Reynolds numbers, heating rate, 
heat input per square foot, total heat input to the vehicle, range, and 
time of flight. The heating analysis considers separately the heating 
of the leading-edge region and the windward side of the vehicle. The 
six deceleration flight paths considered are constant values of altitude, 
angle of attack, G load, heating rate, and free-stream Reynolds number, 
and a case of modulated roll. Limitations on altitude, G load, lift 
coefficient, and Reynolds number have been considered. The following 
conclusions may be drawn: 

1. Of the flight paths considered, all except the constant-heating- 

rate path stayed within the following physical limitations: Altitude 

greater than 0, net acceleration G equal to or less than 10 g r s, and 
lift coefficient less than 0.9 for flight- to circular -velocity ratios 
equal to or less than 2.0 and a vehicle maximum lift-drag ratio of 1.0. 
Entry velocity ratios up to 2.8 are aerodynamically possible. 

2. Changing the type of flight path or the G level of the flight 
path has little effect on the value of the local Reynolds number, al- 
though the value of the free-stream Reynolds number may vary a large 
amount. For the present examples the value of the local Reynolds number 
is about 0.7X10^, a value for which one can be hopeful of achieving a 
laminar boundary layer. 

3 . The flight path with the minimum total heat input is the one with 
the highest allowable G load. 

4. With a laminar boundary layer and a maximum deceleration of 10 
g ! s, entry velocities of twice circular and higher appear aerodynamically 
and thermodynamically feasible for a vehicle with a maximum lift -drag 
ratio of 1.0. About 85 percent of the total heat input to the vehicle is 
to the leading-edge region. The ablated weight is estimated to be about 
5 percent of the vehicle gross weight to decelerate from twice circular 
to zero velocity along a constant 10 G path. A constant-angle -of -attack 
path with an initial G load of 4.25 g’s has a heat input 4.9 times 
higher than that of the constant 10 G path. 

5. For an initial velocity of twice circular, incorporating landing 
point control into a constant 10 G path by controlling the flight path 
in the vicinity of circular velocity increases the heat input about 32 
percent . 


6. For a delta wing vehicle for which the useful volume is the 
volume in the wing, the heat input per square foot varies directly with 
the one-half power of the average fuselage density and inversely with 
the one-half power of the G load and resultant force coefficient. For 
vehicles of this type, with the same sweep back, the total heat input 
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varies inversely with the one-half power of the G load, the one-third 
power of the gross weight, and the one-sixth power of the vehicle average 
density and drag coefficient. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 7, 1961 
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APPENDIX A 
SYMBOLS 


reference area for aerodynamic coefficients, plan area of 
vehicle (see fig. 3), sq ft 

projected frontal area of vehicle (see fig. 3), sq ft 

span of vehicle (see fig. 3), ft 

effective heat-absorbing capacity of ablation heat protection 
system, Btu/lb 

aerodynamic drag coefficient, - — 

I pv% 


aerodynamic drag coefficient at zero lift 
friction coefficient, dimensionless 

component of aerodynamic lift coefficient in vertical direc- 
tion, i 

ipV^A 


total aerodynamic lift coefficient 
pressure coefficient, 2(p 0 - p)/pV 
specific heat at constant pressure, Btu/ ( slug) (°R) 

aerodynamic re suit ant -force coefficient 
2W 

c L 

CdBA 

23W 

root chord of vehicle, ft 

leading-edge axial force coefficient, 

/l 2 

leading-edge axial force/— pV A^ 
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C D 

C L 

C N 

c * 

D 

d 

F 

f 

a 

g 

H 

h 

^w 

h 

J 

«J 

k 

L 


leading-edge drag coefficient, 
leading-edge drag^j- pV^A^. 

leading-edge lift coefficient , 
leading-edge lift pV 2 A f 

leading-edge normal force coefficient, 

/l 2 

leading-edge normal force/— pV A^ 

general wing chord, ft 
aerodynamic drag, lb 

thickness of vehicle (see fig. 3), ft 
see column 19 of table I 

linear distance along streamline from leading-edge stagna- 
tion point to end of leading-edge region (beginning of 
windward side) (see fig. 3), ft 

net acceleration, g ! s 

acceleration due to gravity, essentially the planet surface 
value, ft/sec^ 

see column 19 of table I 

altitude, ft 

total enthalpy at edge of boundary layer, Btu/lb 
wall enthalpy, Btu/lb 

heat-transfer coefficient, Btu/(sq ft)(sec)(°R) 
mechanical equivalent of heat, 778 ft-lb/Btu 
see eq. (7a) 

see columns 11, 17, 18, and 19 of table I 
component of aerodynamic lift in vertical direction 
total aerodynamic lift, lb force 
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l/d 


(L/D) 


max 


l 


M 

m 


n 


P 

Pr 

Q 

% 

% 

q. 

q. 

4 

4 

q l,A=0 


R 


£ 


Re 

Re 


l 


lift -drag ratio 
max 

distance along leading edge of wing (see fig. 3), ft 
Mach number, dimensionless 
entry vehicle gross mass, slugs 

exponent defining variation of altitude density with V for 
general flight path of table I 

static pressure, lb force/sq ft 

Prandtl number, dimensionless 

total heat input, Btu 

total heat input to leading edge, Btu 

total heat input to windward side, Btu 

heat input per square foot, Btu/sq ft 

dimensionless function Q of ref. 4, proportional to total 
heat input per unit area for laminar flow 

dimensionless function q of ref. 4, proportional to heating 
rate for laminar flow 

heating rate, Btu/(sq ft) (sec) 

q^ at 7^=0, maximum heating rate on swept leading edge (see 
fig. 3 for A), Btu/(sq ft) (sec) 

range, ft 

leading-edge radius, ft 

gas constant, 1716 ft -lb/ (slug) (°R) 

free-stream Reynolds number, dimensionless 

Reynolds number at outer edge of boundary layer, dimension- 
less 



radius of planet, ft 
Stanton number, dimensionless 
total temperature, °R 
static temperature, °R 
flight velocity, ft/sec 

component of velocity normal to leading edge (see fig. 3) 

component of velocity normal to any surface (see fig. 3) 

ratio of flight to circular velocity, V/ -/gr, dimensionless 

dimensionless velocity at which change from constant G to 
constant a type of flight path occurs 

entry vehicle gross weight. Earth lb 

wing loading. Earth lb/sq ft 

ablated weight. Earth lb 

distance from leading edge in stream direction, ft 
coordinate systems on vehicle (see fig. 3) 
angle of attack with respect to free stream 
cross-flow angle of attack (see fig. 3) 

exponent describing exponential variation of density with 

altitude in equation = e“^, ft ((3“^ = 2.35x10^ ft 

P SL 

for Earth) 

ratio of specific heats, dimensionless 

angle of surface with respect to free-stream direction 

emissivity, dimensionless 

roll angle (see fig. l) (zero for lift in vertical plane di- 
rected upward) 

conductivity, Btu/(ft) ( sec) (°R) 
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A geometric leading-edge sweep angle (see fig. 3) 

A c effective sweep angle (see fig. 3) 

A leading-edge radius coordinate angle (see fig. 3) 

[i viscosity,, slugs/(ft) ( sec) or (lb force) (sec)/sq ft 

v exponent defining variation of laminar heating rate with wing 

sweepback angle 

p atmospheric density, slugs/cu ft 

average entry vehicle density, Earth lb/cu ft 
p^, average entry vehicle density, slugs/cu ft 

a Stefan-Boltzmann constant, 0.48X10“-^ Btu/(sq ft) (sec) (°R^) 

t time, sec 

cp flight-path angle with respect to local horizontal 

(see fig. l) 

cd leading-edge surface angle (see fig. 3) 

Sub scripts 

aw adiabatic wall 

e at edge of boundary layer 

i beginning of deceleration flight path 

l leading edge 

max maximum 

min minimum 

n general subscript for summation in appendix B 

SL sea level 

s stagnation point 

^Plight quantities of M, p, it, Re, t, t/T, V, y, p, and p with 
no subscript or superscript are free-stream values. 



w 


windward side of vehicle 


w wall 

v vertical component of flight velocity 

1,2 stations along deceleration flight path 

Superscript i 1 

* reference condition in "reference temperature" or "reference 

enthalpy" method of calculating heat transfer 


flight quantities of M, p, (X, Re, t, t/T, V, y 9 p, and p with 
no subscript or superscript are free-stream values- 
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VEHICLE AERODYNAMICS 


The vehicle aerodynamics are estimated "by assuming hypersonic flow 
and using the approximation to the leading edge of the vehicle shown in 
figure 3. 


The pressure coefficient on windward surfaces is in general 

2 

c p = (r + i)sin 2 s = (r + 1 ^\r) 


(Bl) 


and for a surface on the leading edge 

Cp = (r + l 00 ^ 00 ■ ct c )|cos(co - <x c ) (B2) 

In equation (B2) if Cp is positive, the surface is windward, if 

negative, the surface is leeward. The pressure coefficient on leeward 
surfaces is assumed to be zero (rather than the value given by eq. (B2)). 
From figure 3 



sin 2 a + cos 2 a cos 2 A = 1 - cos 2 a sin 2 A 


and 


a 


c 


arc tan 


tan q 
cos A 


(B3) 


(B4) 


Equation (B2) hence becomes 


Cp=(y+l)(l- cos 2 a sin 2 A) cos^co- 


arc tan 


tan a 
cos A 


COSTCO - 


arc tan 


tan ocA 
cos A/ 


(B5) 


For a leading edge composed of multiple flat surfaces of equal width, the 
axial force coefficient of the leading edge using the projected frontal 
area of the vehicle as a reference area is 
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i 


I> P,n 


COS CD. 


n 


C A = 


n 

2 cos ^ 


Similarly the normal force coefficient for the nose is 


(B6) 


n 

£ °P,n sin % 

1 

cos A 

1 


n 

2 COS “n 


(B7) 


The contribution of the 
suits of equations (B6) 


leading edge to the drag and lift, from the re- 
and (B7 ) , is 

cos a + c^ sin a (B8) 

C L = C N cos a " c a s ^- n a (B9) 


The total lift and drag coefficients of the vehicle based on the refer- 
ence area of the vehicle (see fig. 3) are then 


2 

•X,a = (T + l)sin a cos a + 


(BIO) 


A f 

C D = (r + l)sin 3 a + c D -j- 


(Bll) 


where the first terms on the right side of the equations are the contri- 
bution of the reference area A to lift and drag coefficients, and the 
second terms are the contribution of the leading edge. 

If one considers a family of vehicles having the general plan form 

shown in figure 3, that is, having the same sweepback angle but various 

values of c/d, or what is equivalent, various values of (L/D) , then 

1 mux 

the following relation may be derived for low angles of attack: 

A f 

C D ~ C D,o = c D,o T 


TnnT- 
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where Cj,^ 0 is a constant for a given sweepback angle. From equation 
(El) 



and 


so that 


A f = bd 


c 

C-p — = 2 cj^ q = Constant 


(B12) 
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APPENDIX C 


WINDWARD-SIDE HEATING ANALYSIS 

A simple analytic expression for the laminar heat transfer to the 
windward side of the vehicle at high velocities is developed herein. 

Off the stagnation point for moderate angles of attack the heat 
transfer to the windward side is analyzed using an approximation to the 
"reference temperature" (ref. 9). The relations required are: 


4 = h ^aw “ 

(Cl) 

St - h 

p*V c 

H e p 

(C2) 

Q 

St = (Pr*)“^/^ 

(C3) 

C f - °- 66 ^pW e x 

(C4) 

1 + v^ Z v j - vp? 2 J pJ 

(C5) 

P*C* 
Pr* = — 

K* 

(C6) 


t* = t e + 0.5(t^ - t e ) + 0.22(t aw - t e ) (C7) 

Assuming t^ « t aw and combining equations (Cl) to (C5) yield 


q. = 


0.166 


JPr* 


0.167 


4 


p*Vu> 


V‘ 


(C8) 


The three terms in the numerator of the square root term of equation 
(C8) may be evaluated as follows, as an approximation of Southerland's 
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viscosity equation for high temperatures (ref. 17, eq_. (A3)): 

,3/2 


u* = 2.27xl0“ 8 -4£ — —— 2.27X10" 8 t* 1 / 2 

t* + 198 


(C9) 


The reference density is given by 


p* = p fk (C10 ) 

At moderate angles of attack, by Newtonian aerodynamics, 

v e = V cos a (Cll) 

The terms t* and P g /p in equations (C9) and (CIO) must now be 
evaluated. In equation (C7) for low and moderate angles of attack, t e 
and (t^ - t e ) are small compared with (t aw - t e ) and t aw , so that equa- 
tion (C7), combined with equation (C5), becomes 


t ^ 0 22 t - 0 - U V Pr * V 

t 13 aw " C p J 


(C7a) 


From hypersonic aerodynamics 

/P e 

«P,e - ( T 


■ 1 )t^ 


(y + l)sin 2 cx 


(C13) 


where at super circular speed (M > 25) and for moderate angles of attack 
P e /p » 1.0. At the low angles of attack where the assumption produces 

maximum error, the heat input to the windward side of the vehicle is gen- 
erally small compared with the heat input to the leading edge (see fig. 
14:) . The error in estimating the total heat input to the vehicle can 
thus still be small at small angles of attack. Then, since 


M 


2 



equation (C13) becomes 


£e _ x + 1 V 2 

p ~ 2 £t 


. 2 

sin a 


y + 1 r 

2 ftt 


. 9 

sm a 


(C14) 



Making the appropriate substitution in equation (C8) gives 


q. = 


0.166 


y 2.27X10' 


j0 . 75 pr *0 . 292 ”1 / lo. 11 1 2 


* 


+ 1 pV 2 cos q sin 2 a tt2 


V gr (C15) 


The following values are assumed to be constants: 


Pr* =0.71 

T = 1.2 

ft = 1716 ft -lb/ ( slug) (°R) 

C P = l) = 13 - 2 ®tu/ (slug) (°R) 

J = 778 ft-lb/Btu 

From the hypersonic aerodynamics of the configuration, equation 
(Bll), and equation (6) 


A cos 9 


qr = — pV^ sin^a cos a 
2k c 


(C16) 


where 


k c - 1 + 


AfCp 


V-l/2 


A(x + l)sin^a cos a/ 
so that equation (C15) becomes, with the assigned constants, 


(C17) 




q Tr = 1.5X10 k gr* WA , . , 

w c e TAx cos 0 1 


■{* 


V 2 


(C 18 ) 


The absolute value signs are used to indicate that the heating is real 
and positive irrespective of the sign of L and cos 0. By use of the 

value of — given in column 6 of table I, one may solve for L in 
g ClT 

equation (4b) and substitute in equation (C18) to arrive at the results 
given in column 10 of table I. For example, for the constant Re case 
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with 0=0 



(C19) 


2 

The term k c is the ratio of the lift produced by the vehicle wind- 
ward side to the sum of the windward-side and leading-edge lift . Approx- 
imately, for angles of attack from 0° to 15° 


d( sin a) 


where dc^/dCsin cc) is approximately constant and depends on the sweep 

angle. For the present configuration with A = 60° and 
dc^/d^sin a) = 0.88, A^/A = 0.569, a = 5.3°, and k c = 0.535. 


The term k c , in general, varies over the flight path. For super- 
circular entry velocities a conservative estimate of the heat input is 
obtained by evaluating k c at the maximum velocity. This is done in 

the present calculations. The value of k c approaches unity and con- 
stancy for small values of A^/A or high values of (L/D) max . 
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APPENDIX D 


METHOD USED FOR CALCULATING RATIO OF LOCAL 
TO FREE-STREAM REYNOLDS NUMBER 

The procedure for calculating the local Reynolds number at super- 
sonic speeds for cones and wedges with unswept edges is presented in 
reference 18. A relation for calculating the local Reynolds number at 
hypersonic speeds for the case of a flat plate at angle of attack with a 
swept blunt leading edge is developed herein. The nomenclature is illus- 
trated in figures 3 and 18. 

At a station sufficiently far downstream from the leading edge, 
where the surface static pressure for a blunt body is approximately that 
for a nonblunt body, the ratio of local to free-stream Reynolds number 
per foot or for the same length may be written 

^ Is. JL (Dl) 

Re p V p e v ; 

In equation (Dl) the velocity ratio is taken as the Newtonian value. 


= cos a 


(C9), 


The viscosity approximation is the same as that used in equation 


u = 2.27X10 -8 t 1 / 2 
^e e 


where for this case 


= t s - 


where 


t s = t + 


cos 8 A c cos^A c 
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Hence equation (D3) becomes 


= 2.27X10 


-8 


r-i 

V cos A_/p.\ 2 


V2JC 


P 


The term P e /p in equation (Dl) may he written 

Pe _ Ps Pe 

~ “ ~ pI 


and equation (Dl) then becomes 
Re 


/ P s 

P cos aV2JC p — 

Re 2.27X10" 8 V-/ir cos A \P 


The ratio P e /p s may be written by use of the perfect gas law as 


VP / 

C v S' 


5-r 

2 


Haf-eri) 


i/r /D .x -i/r 


Hypersonic aerodynamics gives (see also eq. (C14)) 


(D6) 


(D7) 


(D8) 


(D9) 


= 1 + 


r + l Yl 

2 &t 


sin a 


(DIO) 


and 


h = ! + v£ cos2 » nil v£ oos2 

p 2 &t G 2 fct c 


(Dll) 


where the approximation is valid for sweep angles less than about 87°. 
From figure 3 (see also eq. (B3)) 


eos^A c = (l - sin^A cos^a) 


(D12) 


The density ratio across the shock at the cross-flow stagnation point is 

P, 


l s _ x + 1 

p ” r - i 


(D13) 
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With the appropriate substitution, equation (D8) gives finally 



(x + l\ 

2Rt , . ? 

p cos a-/2JC p ^ r -J 

(r + l)v 2 gr 

2.27X10“® V-y/gr(l - sin 2 A cos 2 ^ 1 ^ 2 

o 2 

1 - sin A cos a 


5-r 

2r 


(D14) 


The following values were assumed, in the present calculations : 

T = 1.2 

ft = 1716 ft -lb/ (slug) (°R) 

C P = i) ' “ 13 • 2 ®tu/ ( slug ) ( °R ) 

t = 500° R 

(j. = 0.37X10“® (lb force) (sec)/ sq ft 
J = 778 ft-lb/Btu 

The term V is the independent variable, and a is determined by the 
lift coefficient (column 7 of table i) and the vehicle aerodynamics. 

This calculation takes no detailed account of the boundary-layer 
history on the leading edge of the vehicle. 

Typical curves of the ratio of local to free-stream Reynolds number 
are shown in figure 18 as a function of angle of attack for several 
values of V. Both V and angle of attack have a strong effect on the 
Reynolds number ratio. 


fc= 

i 


-mm - ' 
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APPENDIX E 


INTEGRATION OF HEAT INPUT OVER VEHICLE 

The vehicle geometry over which the integration of the local heat 
per square foot q. is carried out is shown in figure 3 together with 
the symbols and coordinate systems used. The constant Re case is used 
as an example . 


Vehicle Geometry 

The following approximate relations can be derived from the vehicle 
geometry : 


. 1 . _ W m 

2 tan A W JK in/A 


A r = bd = — . 
1 p f \tan Aj 



m 


m 


1/2 


_f 

A 


pj. \tan AJ 


2d = 2(W/A) 5 / 2 

c Pf (W tan A) 1 / 2 


(El) 



b = 22 cos A 


c = 




tan A \ 1//2 

w~) 


W 

% A 


— 2p . 


or 


— = 2p« 

%A P f 
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Integration Over Leading-Edge Region 


The integration over the vehicle is divided into two parts, that 
over leading-edge regions, and that over the windward side of the ve- 
hicle. To illustrate the integration over the leading-edge regions, q. z 

may he taken from column 15 of table I (hence the subscript 15 on k) or 
the integration of equation (37) and written as 


^7 


= k-^g £l - sin 2 A cos 2 oj ^ cos A 


(E2) 


where for constant Re 


k i5 - 12 - 1 x1 ° 3 -4/2a c 1 


vf - 5 - v |- 5 

2T5 


N D _ 

^ CR 


via 


The general form of the integration desired is (see fig. 3) 

•7 r it/2 


Q z = 2 



R*r dA d7 


0 •4.(jt/2-a c ) 


^7% 


(E3) 


Carrying out the interior integration and using equations (El) and (E2) 
yield 


/ 


'/ 2 


Q z = k 15 d | [l - sin 2 A cos 2 a] (l + 


cos a c )d7 


where from figure 3 


, tan a 

= arc tan r 

c cos A 


Equation (E4) may be written 


/ 


Q z = 2k-|_5d J k z cos v A d7 


(E4) 


(E5) 


(E6) 



E-1001 


45 


where 


k 2 = 




sin^A costx 


& 


1 + cos arc tan 


/tan ccA ) 
\cos A // 


2 cos v A 


(E7) 


For a = 0, k z = 1.0; for a = 12°, v = 1.5, and A = 60°, k z = 1.05; 

thus, for the integrations over the flight paths for which the angle of 
attack is small, k z is assumed to he unity. For the constant a case, 

which is used at higher a, k z is a constant and may he evaluated. With 
k z a constant and hy use of equation (El) , equation (E6) can he written 


Qt - 


k 15 k 2 d 


COS 


(v-1)a 


/•b/2 
/ dh 

Jo 


(E8) 


so that finally integration gives 

Q z = i k^cjk^bd cos^-^A = i k^k^A^ cos^ v “-*-^A (E9) 

Writing this result in terms of the gross weight and wing loading, hy 
using equation (El) , yields 


0 1_ _ k 15 k Z / w/A W 2 

W " \W tan A/ 


(E10) 


or 


% = k 15 k 2 c^ (v ~ 1)a / m /A N 1 / 2 
m p ^ yn tan A / 

and substituting for k-^g (eq. (E2)) and for W/RjjA from equation (El) 
results in 
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or 


% 

~ ” ' -9 3/2! (v-l)-v 

21.9x10 mgr 7 k z cos A 


= vs ( , m / A .A 

yp m tan Ay 


1/2 


V 2,5 _ v 2 ' 5 
1 2 


2.5 


\Pf m tan A ) U D /TgT 

Cd ^Vk 7 ) 


"R 


Comparing equations (E2) and (Ell) yields 

Q 


w 


V 2 k z cos^ v_ ^A 


(l- sin^A cos^a)^ 2 cos A /Pf W ' t:an ^ ^ 

W/A 



(E12) 


% A 


Equation (E12) results from an integration over the vehicle and hence 
is true for all flight paths. It is listed in column 17 of table I and 
a similar relation is given in table II. 


Integration Over Windward Side of Airplane 

For this case q w may be taken from column 16 of table I (hence 
the subscript 16 on k) or the integration of equation (38) and may be 
written as 


% = 



(E13) 


where , 


for example, for the constant 
k 16 =8.2X10 3 k c ^| 


Re case for V > 1, 

- 1 1 %e> 5/ f 1 

k 



The integration desired is then (see fig. 3) 


^w = 



d| dq 

V | - t] tan A 


(E14) 


E-1001 


47 




Comparing equations (E13) and (E19) gives 


8 



(E19) 


(E20) 


Equation (E20) holds for all flight paths and is given in column 18 of 
table I. A similar result is given in table II. 
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From section AA of figure 3 the ratio f/c may be written approxi- 
mately for small angles of attack as 


f n i _ 4 *] 

c 4 c [ cos A ~ ir ac J 


and. by use of equation (E5) 


(E21) 


f 

c 


A 

4 




1 

cos A 


- — arc tan 
n 


tan a 
cos A 


(E22) 


For the configuration used as an example in this report to generate 
numerical values, % = 2.06, c - 14.5, and A = 60°. For a = 0° and 

5.3°, f/c is 0.34 and 0.31, respectively. The corresponding values of 
(eq. (E17)) are 0.48 and 0.49, respectively. Note that k w does not 
vary significantly for small changes in a, and hence it is assumed to 
be constant. A value of 0.49 was used for the calculations presented. 

The most desirable form of some of the terms in equations (Ell) and 
(E19) depends on what information is known. The following useful alter- 
nate forms may be derived if the parameter W/A is eliminated by use of 
one of the equations (El) : 



The term c/d depends on the vehicle sweep angle and maximum lift-drag 
ratio. For the present examples with A = 60°, c/d = 3.51 and 21.7 for 
maximum lift -drag ratios of 1 and 2, respectively. Still another inter- 
esting form of the results of columns 17 and 18 may be written using the 
results of equations (E23), (E24), and (B13) . For vehicle operation at 
low angles of attack where ~ ~ , which is approximately the 

same assumption made to arrive at the results of columns 17 and 18, 

VS f^w/T A ) ' C J /6 W l/3 Pf /6 ( 2c D A > 1/3 


(E25) 
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0 W tan a) 1/4 = cj/M/^pl/e tani A a_ (E26) 

where is constant for a given sweepback angle. 

An example of the significance of this result is the following. For 
a family of vehicles having the same sweepback angle , gross weight, and 
fuselage density, the total heat input to the vehicle leading edge varies 
as 


% 


i 

pTe 


(E27) 


(by the results of column 17 and eq. (E25) for kj = Constant) . On the 

other hand (by column 15) the heat input per square foot to the leading 
edge varies as 


~ -J/z (E28 > 

C D 

Thus the change in the total heat input to the vehicle is much less 
sensitive to changes in Cp (or (V D ) max ) than is the heat input per 
square foot. 
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APPENDIX F 


CASE OF CONSTANT ANGLE OF ATTACK 


The determination of the term 


V dcp 
g 


is more difficult for the case 


of constant angle of attack than for the other flight paths. However, 
the solution of this case permits a check on the accuracy of the results 
of the present analysis by comparison with the numerical integrations of 
references 4 and 13. The equations to be manipulated are repeated here 


1 dV 'J. 
g dx ' U CL 


G-^=1-V 2 +-j& 
Cr g dx 

C R pV 2 Agr 


C PS VA dp 


2(3W d y 


Using the mathematical identity 


V d£ _ V d£ dV 
g dx g d y dT 

and equation (3) permits equation (4a) to be written 


1 - V 2 


^ 1 + 


V d£ 


Solving equation (5) for p, by use of (F2) for G, yields 


P = C p(p 


where 



Differentiating equation (19) gives 



where 

C D gA 

C 9“W 


(F4) 


(F5) 


(F6) 


The term dcp/dV can be eliminated between equations (F3) and (F5) to 
yield a single differential equation for which a solution must be found. 
Herein, however, a solution to equations (F3) and (F5) was found by 
successive approximations. For a first approximation the flight-path 


angle is assumed to be a constant 



Then from equation (F3) 


and 





and equation (F5) becomes 


dqp _ 4C cp C p _ 4 

dV ” V 3 " prv 3 jjj 


(F7) 


(F8) 


(F9) 


The value for dcp/dV given by equation (F9) may now be introduced into 
equation (F3) to obtain a second approximation. Thus, 


P = 



(F10) 
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where 


3r(L/D) 2 


Repeating the steps of equations (F7) to (F9) results in 



2k A 

v 5 v ? 2 



(in) 


(F12) 



Higher order approximations may "be obtained. If this degree of approxi 
mation is accepted, equation (F2) for G becomes 


where 



(F15) 


(F16) 


With G known, p can be determined from equation (5), and this is the 
result given in table I. 


From equations (19) and (F12) 


where 


2H 

(3rV2 k 



k(l - V 2 )1 

V 2 + k J 


(FI 7 ) 


(F18) 
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From equation (Fl), by use of equations (3) and (F14) , 


= _ 1^1 (i 

g dT F 



which is the result given in 


column 6 


of table I . 


(F19) 


I 



APPENDIX G 


COMPARISON WITH OTHER ANALYSES 

The results of reference 4 are essentially for a three-dimensional 
stagnation point, or j = 1 in equation (7a). The result of column 10 
of table I for constant angle of attack for a three-dimensional stagna- 
tion point can be written in the modified form 


4= V? 14.9 

Vl7d V f 


m 


(Gl) 


From equation (36) of reference 4 for Earth, and in the nomenclature of 
this report. 


q = — r 

590j- 

V c 


m 


c d r n a 

Substituting equation (Gl) in (G2) gives 


(G2) 


q = 0.202 


Vl - V 2 V 2 




(03) 


The equation from column 15 of table I for a three-dimensional nose is, 
in modified form. 


q = -y/2 12.05Xl0 3 Vg 



m_ 

C D% A 


arc sin. V 




■1V n 


(G4) 


From equation (39a) of reference 4, and in the nomenclature of this 
report , 


q. = 


15,900 




(G5) 


m 


CdR n a 
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Evaluating (G4) for V-j_ = 0.99 and Vg = 0, the values used in reference 
4, and substituting in equation (G5) yield 


i = 3.91^| 


(G6) 


As a point of interest, the total heat input per square foot q or q 
is quite sensitive to the initial velocity for this case. For example, 
evaluating equation (G4) from V-^ = 1.0 to Vg =■ 0 gives for equation 

(G6) 


q. = 


4.76 -^”1 


(G7) 


which is an increase of 22 percent over the case for Vj_ = 0.99. In the 

present calculation the term F is evaluated at the highest value of 
V . If the flight path_is divided into several parts, F is evaluated at 
the highest value of V in each part . 

Reference 4 also presents the horizontal component of the G load, 
Gjj. For small flight-path angles, an assumption already made in the 

present analysis, 


C-n 

Gjj ~ G J7— (G8) 

°R 

With the result given for G in column 4 of table I for the constant a 
case. 



(G9) 


Part of the difference in heating rate between the present analysis 
and reference 4 is due to the use of a slightly different value of the 
constant in equation (7a) . 
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Point functions j 

Any planet 


Plight path 

Constant of 
flight path 

Altitude density, 
slugs/ cu ft 

0 load, 

g's 

Plight-path angle, 
v, 

radians 

IS* 

Lift coef- | 

flclent. 

Free- 

stream 

Reynolds 

number, 

uRe/l 


Modulated angle of attack; 0-0 
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Auxiliary relations or definitions 



Two-dimensional leading edge. 

Where two equations are entered at one location, the applicable equation is the one that yields real numbers. 



FLIGHT-PATH CHARACTERISTICS 
are always positive. ] 


Functions Integrated over flight path and over vehicle 


Heat input to leadlng-edge c region for Heat input to windward side for 

a ”> 0° to 30 6 , Qj/W, a - 0° to 30° , Q^/w, 


Btu/lb of vehicle gross weight 
(e) 


Btu/lb of vehicle gross weight 

(f) 


Modulated angle of attack;; 0 = 0 
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Auxiliary relations or definitions 


{l- sin2 A cos 2 a} /2 Jl + cos[arctan(g|^)]j k w - (l + |) 7 - (|) 7 ~ |(|) 7 


7 = & tan 


wr~ ) = (cr tan A j w p; 


(see also eq. (E25)) 


l, is given by column 15. 


} W tan A V* . tan^A „l/5 p l/6 

(I) 

(see also eq, (E26) ) 
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(a) Lift coefficient. 


Figure 2. 


- Aerodynamic characteristics 
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(b) Drag coefficient. 

Figure 2. - Continued. Aerodynamic characteristics. 
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(f) Lift -drag ratio. 

Figure 2. - Concluded. Aerodynamic characteristics. 
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Nondimens ional total heat Nondimens ional heating rate at Horizontal component of 
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Present analysis 


Ref. 4, fig. 10 
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Present analysis 


Ref. 13, fig. 4 
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Figure 6. - Constant-angle-of-attack flight path for Earth atmos- 
phere. Angle of attack, 5.3°; initial net acceleration, 10 g T 3; 
initial velocity ratio, 2.0; maximum lift-drag ratio, 1.0; entry 
vehicle gross weight, 10,000 pounds; average entry vehicle den- 
sity, 20 pounds per cubic foot; geometric leading-edge sweep 
angle, 60°. 
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Figure 7. - Constant-altitude flight path for Earth atmosphere. 
Altitude, 162,000 feet; initial net acceleration, 10 g's; 
initial velocity ratio, 2.0; maximum lift-drag ratio, 1.0; 
entry vehicle gross weight, 10,000 pounds; average entry vehicle 
density, 20 pounds per cubic foot; geometric leading-edge sweep 
angle, 60°. 
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Ratio of flight to circular velocity, V 

Figure 7. - Concluded. Constant-altitude flight path for Earth 
atmosphere. Altitude, 162,000 feet; initial net acceleration, 
10 g’s; initial velocity ratio, 2.0; maximum lift-drag ratio, 
1.0; entry vehicle gross weight, 10,000 pounds; average entry 
vehicle density, 20 pounds per cubic foot; geometric leading- 
edge sweep angle, 60°. 
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Figure Q. - Constant free- stream Reynolds number flight path for Earth atmos 
phere. Free -stream Reynolds number, 7.3X10^; initial net acceleration, 

10 g’s; initial velocity ratio, 2.0; maximum lift-drag ratio, 1.0; entry 
vehicle gross weight, 10,000 pounds; average entry vehicle density, 20 
pounds per cubic foot; geometric leading-edge sweep angle, 60°. 
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igure 8. - Concluded. Constant free-stream Reynolds number flight path for 
Earth atmosphere. Free-stream Reynolds number, 7.3X10 6 ; initial net accel 
eration, 10 g’s; initial velocity ratio, 2.0; maximum lift-drag ratio, 1.0 
entry vehicle gross weight, 10,000 pounds; average entry vehicle density, 
20 pounds per cubic foot; geometric leading-edge sweep angle, 60°. 
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Figure 9. - Concluded. Constant-net-acceleration flight path for 
Earth atmosphere. Net acceleration, 10 g f s; initial velocity 
ratio, 2.0; maximum lift-drag ratio, 1.0; entry vehicle gross 
weight, 10,000 pounds; average entry vehicle density, 20 pounds 
per cubic foot; geometric leading-edge sweep angle, 60°. 
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Figure 10. - Concluded. Constant-heating-rate flight path for 
Earth atmosphere. Maximum leading-edge heating rate, 900 Btu 
per square foot per second; initial net acceleration, 10 g*s; 
initial velocity ratio, 2.0; maximum lift-drag ratio, 1.0; 
entry vehicle gross weight, 10,000 pounds; average entry 
vehicle density, 20 pounds per cubic foot; geometric leading- 
edge sweep angle, 60°. 
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Figure 11. - Modulated roll flight path for Earth atmosphere. Initial roll 
angle, 0; initial net acceleration, 10 g T s; initial velocity ratio, 2.0; 
maximum lift -drag ratio, 1.0; entry vehicle gross weight, 10,000 pounds; 
average entry vehicle density, 20 pounds per cubic foot; geometric leading 
edge sweep angle, 60°. 
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Figure 11. - Concluded. Modulated roll flight path for Earth 
atmosphere. Initial roll angle, 0; initial net acceleration, 
10 g's; initial velocity ratio, 2.0; maximum lift -drag ratio, 
1.0; entry vehicle gross weight, 10,000 pounds; average entry 
vehicle density, 20 pounds per cubic foot; geometric leading- 
edge sweep angle, 60°. 
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Figure 12. - Total heat input at Earth for flight paths with j 

net acceleration of 10 g T s at velocity ratio of 2.0. Maxi- 
mum lift-drag ratio, 1.0; entry vehicle gross weight, 10,000 
pounds; average entry vehicle density, 20 pounds per cubic 
foot; geometric leading-edge sweep angle, 60°. 
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Figure 13. - Total heat input at Earth for flight path with initial 
net acceleration of 10 g f s. Maximum lift-drag ratio, 1.0; entry 
vehicle gross weight, 10,000 pounds; average entry vehicle den- 
sity, 20 pounds per cubic foot; geometric leading-edge sweep 
angle, 60°. 
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Figure 14. - Total heat input for constant-angle-of -attack flight path for 
Earth atmosphere. Initial velocity ratio, 2.0; maximum lift-drag ratio, 
1.0; entry vehicle gross weight, 10,000 pounds; average entry vehicle 
density, 20 pounds per cubic foot; geometric leading-edge sweep angle, 
S0°. 
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(a) Modified constant -angle -of -attack flight path. 
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(h) Modified constant -net -acceleration-flight path. 
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(c) Composite flight path. 

Figure 15. - Flight paths with control of landing location. 
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Figure 16. - Total heat input at Earth for modified constant -angle- 
of -attack flight path (fig. 15(a)). Initial net acceleration, 

10 g's; initial velocity ratio, 2.0; maximum lift -drag ratio, 1.0 
entry vehicle gross weight, 10, 000 pounds; average entry vehicle 
density, 20; geometric leading-edge sweep angle, 60 . 
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Figure 17. - Total heat input at Earth for modified 
constant-net-acceleration flight path (fig. 15(b)) 
and infinite range. Initial net acceleration, 10 g 1 £ 
initial velocity ratio, 2.0; maximum lift-drag ratio, 
1.0; entry vehicle gross weight, 10,000 pounds; aver- 
age entry vehicle density, 20; geometric leading-edge 
sweep angle, 60°. 
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Figure 18. - Ratio of local to free-stream Reynolds number for 
swept blunt leading edge for Earth atmosphere. Leading-edge 
sweep angle, 60; specific-heat ratio, 1.2. 
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